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^ . Abstract 
X 

I ' Inclusive event shape variables have been measured in the Breit Frame for neutral 

current deep-inelastic positron-proton scattering using the HI and ZEUS detectors at 
r~| ! HERA. The variables thrust, jet broadening, C-parameter, jet mass and two kinds of 

differential two-jet rate have been studied in the kinematic range 7 < Q < 100 GeV. The 
k> ' ^ dependence of the shape variables have been compared with QCD applying power cor- 

, rections proportional to to account for hadronisation effects. The concept of power 

I corrections is tested by fitting the strong coupling constant Og and a non-perturbative 

parameter cSp-i. 



1. Introduction 



A recent revival of interest in event shapes has been prompted by theoretical developments in the 
understanding of hadronisation [Q. These approaches extend perturbative QCD calculations 
into the region of low momentum transfers using approximations to higher-order graphs; the 
hadronisation contributions are expected to go as where p is an integer power (1 or 2). 

The scale at HERA can be varied over four orders of magnitude allowing these corrections to 
be studied in detail. The data are compared to theoretical expectations for power corrections, 
characterised by as{Mz) and an effective coupling ap_i(/i/) specified at the infrared matching 
scale /X/ ~ 2 GeV , and Next-to-Leading Order (NLO) perturbative QCD (pQCD) calculations, 
determined by as{Mz). 

A suitable frame of reference in which to study event shapes is the Breit frame since this 
frame maximises the separation between the current jet from the struck quark and the proton 
remnant. Viewed in this frame the exchanged gauge boson is purely spacelike and the incoming 
QPM quark is back-scattered with equal and opposite momentum {Q/2). The Breit frame also 
permits the event shapes from the spacelike process in DIS to be compared directly to the event 
shapes in the timelike e~^e~ process. Two classes of event shapes are studied. 

The first class are those confined to the current region of the Breit frame {z < 0), namely, 
thrust, both with respect to the Breit frame axis (r^) and the thrust axis (Tmax), jet broadening 
with respect to the Breit frame axis (B), jet-mass (p), and the C-parameter (C). Thrust (r^) 
and jet broadening (B) are measured with respect to the Breit frame axis and are therefore 
specific to DIS; the remainder can also be defined in e~^e~ experiments. At the hadron level, the 
HI result 0] assumes the hadrons have their true mass, the ZEUS result assumes the hadrons 
have zero mass; this difference in definition only affects the jet-mass which can be seen in Fig. 
Here the data are compared and the other variables show agreement between the two data sets. 
ZEUS bin in kinematic ranges of both x and Q, so in some regions ZEUS have two data points 
measured over the same range in Q but a different range in x; HI integrate over x. In order to 
keep the variables infrared safe, the total energy in the current region has to exceed 10% (3%) 
of the total available energy, for HI (ZEUS). The exact value of this cut is not critical. The 
second class of variables are the two-jet rates. These are measured over the full phase space 
of the Breit frame in which (n + 1) jets are searched for; the +1 denotes the proton remnant. 
Both the kt (or Durham) and the factorisable JADE scheme are used in which a jet merging 
criteria based on distance measures yij between pairs of particles is defined. The event shape 
y denotes the value of yij at which the transition (2 + 1) ^ (1 + 1) jets occurs. 

The mean values of these event shapes can be interpreted theoretically as 



(F) = {Ff"'' + (1) 

where 

■c^m (2) 



{Fy^'' = c,,f{x, Q) ■ as{Q) + C2,^(x, Q) ■ 



with ci,i?(x, Q) and C2,f{x,Q) being calculable coefficients of pQCD, here determined by the 
DISENT NLO program with renormalisation/factorisation scales set at fiR = fip = Q. A 
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Figure 1. Comparison of HI and ZEUS data. Note the x-dependence in the ZEUS data. 
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simple form for (_F)p°™ oi A/Q oi B /Q'^ has been ruled out by HI data 0. A more sopisticated 
approach introduces 



pow 



ap_i(/i/) - - ( In — + 0.45 + - ) al{ii 



(3) 



where p is the power, ap is a calculable F-dependent coefficient and A^' = 2M./ti ~ 0.95 |Q is 
a two-loop level refinement to the calculations called the Milan factor. 



2. Results 



Fits of Eq dlD to ZEUS data are plotted in Fig. |^. The shaded band indicates the theoretical 
uncertainty on (F)^*^"^* from varying the renormalisation scale by a factor of 2 to = 2-Q and 
/iR = 1/2 ■ Q; this is the dominant uncertainty on the measurement. The fitted values of the 
two parameters as{Mz) and ao{fii) from HI and ZEUS data are plotted in Figs. 3 and 4 as 
contours (theoretical errors such as the scale dependence are not included in the contours). The 
HI result shows the = Xmin + 1 = Xmin + ^ contours for statistical and experimental 

systematic errors combined in quadrature. The smaller contours in the ZEUS result are a 
statistical only fit, the larger ones are the = Xmin + 4 contours for the statistical and 
systematic errors combined in quadrature. Large correlations between as{Mz) and ao{fii) are 
seen for (r^) but less for the other variables. Both experiments fit an approximately universal 
ao{fj,j) ~ 0.5 ± 0.1%. The fitted values of as{Mz) for the e+e^-type variables (C-parameter, 
jet-mass and thrust w.r.t. thrust axis) are compatible for the ZEUS result and remarkably close 
for the HI result. Including the DIS specific variables (thrust and broadening w.r.t. Breit axis) 
sees a larger spread in the fitted values of as{Mz)- Possible explanations are missing higher 
order QCD corrections or incomplete knowledge of the power correction coefficients. Fits to 
HI two-jet rate data are plotted in Fig. ^ [Q. The hadronisation correction is seen to be small 
for these variables. For the (y) from the factorisable JADE algorithm, table |l] shows that the 
fitted value of aoif^i) is unreasonably low with the conjectured ap factor of 1.0. Instead, a 
much more reasonable fit can be obtained by allowing the power corrections to become small 
and negative, with an ap factor of —0.2. 



Table 1 

Fits of {y) to HI data showing preference for negative power corrections. 



(F) 


ap 


P 


ttp-i(/i/ = 2 GeV) 


asiMz) 


xVdof 


iVfj) 


1 


1 


0.28 


n +0.006 
U.iUd _o.006 


0.8 


iVfj) 


-0.2* 


1 


n 07 +U.2U 
u-J' -0.21 


nfi +"-'^"^ 

U.iiO _0.009 


0.6 



The power fit for the {y) from the kt algorithm is of special interest because the power p is 
expected to be 2 instead of 1 as it is for all the other variables. The HI data are inconsistent 
with a fit using p = 1. The data are consistent with p = 2 but a more qualitative statement 
cannot be made given current experimental precision on this variable and the fact that the 
theoretical ap factor for this variable is unknown. 
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Figure 2. Power correction fit to ZEUS data. The dashed hne is the NLO (F)'^'^^ contribution, 
with the shaded band showing the effect of varying the renormahsation scale by a factor 2. The 
dotted line is the power contribution and the full line is the sum of the two. The region 

below 9 GeV is excluded from the fit. 
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Figure 3. HI contours showing fitted values of as{Mz) and ao{fii). The plot sliows tlie statis- 
tical © systematic contours for two confidence levels for each of the variables. 



5 



0.7 



> 
(D 
O 

CM 



III 



0.65 - 



0.6 



0.55 - 



0.5 



0.45 



0.4 



0.35 



0.5 



0.09 



ZEUS Preliminary 1995-97 



- 






• btat — only tit, U>y beV 
A Stat + Sys tit, Q>9 GeV 
(x' = Xmin + 4 contours) 


















- 














Po 

^ c 








\ ^f\\ 


World Average 

1 






\ ^ \ 

\ 1 \ 1 



0.1 



0.1 1 



0.12 



0,13 



0.14 



0.15 
«s(Mz) 



Figure 4. ZEUS contours showing fitted values of as{Mz) and Q;o(/i/). Tlie plot shows both 
the statistical only contours and the statistical © systematic contours for each of the variables. 
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Figure 5. Fits of (y) HI data, for kt (left) and factorisable JADE (right) algorithms. The 
dotted line is the perturbative QCD prediction only, the solid line is this prediction plus the 
power correction. 
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ZEUS have also studied the x-dependence of the power correction fits. Figure || shows the result 
of omitting the lower x bins from the ZEUS data given in Fig. ^ plotted are the statistical 
contours only. Most of the variables display little sensitivity to x. However, a large effect is 
seen in the jet broadening where the fitted values of as{Mz) and ao{fj,i) shift considerably. This 
confirms recent theoretical comment that the ap factor for jet broadening should contain 
x-dependent terms that yet, unknown. 

In conclusion, this new theoretical approach to hadronisation works very well given its simplic- 
ity, although there remain some theoretical questions particularly in the x-dependence of the 
DIS specific variables. 
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Figure 6. Fit to ZEUS data showing a clear shift in the contour for jet broadening due to 
the x-dependence. Plotted are the statistical only contours. 
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